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INTRODUCTION 


Research  under  contract  F4920-79-26  began  in  1978  with  the  intent  of 
gaining  a  better  understanding  of  the  aerodynamic  behavior  of  unsteady  flows. 
The  unsteady  excitation  was  generated  by  a  rotating  elliptical  cylinder 
located  below  and  behind  the  airfoil.  This  produced  an  unsteady  flow  of 
essentially  constant  phase.  Since  1978,  testing  has  been  performed  with 
reduced  frequency,  based  on  half  chord,  varying  from  0  (steady  flow)  to  6.4. 
The  Reynolds  number  varied  from  375,000  to  1,400,000,  but  most  of  the  testing 
was  conducted  at  a  Reynolds  number  of  700,000.  The  metric  surface  used  was  a 
NACA  0012  airfoil,  of  20  inch  span  and  chord,  which  was  set  between  two 
sidewalls  to  ensure  two-dimensionality.  The  details  of  instrumentation  are 
shown  in  figures  1  and  2  and  described  in  depth  in  references  1,  2  and  3.  The 
unsteady  data  taken  during  a  test  is  separated  into  three  parts:  a  time 
invariant  mean  part;  a  periodic  part  of  zero  mean,  but  not  necessarily 
sinusodial;  and  the  remainder  which  has  zero  mean  and  is  not  periodic. 
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The  time  invariant  part,  or  time  mean  part,  f,  is  the  time  average  over  a 
cycle.  The  periodic  part,  f,  is  the  ensemble  average  of  the  cycles  less  the 
time  mean.  The  remainder,  f',  is  the  difference  between  the  signal  and  the 
ensemble  average,  and  is  assumed  turbulent. 

Each  year  since  work  began  in  1978,  a  number  of  significant 
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accomplishments  have  been  made.  Work  began  in  1978  with  studies  of  turbulent 
boundary  layers  on  a  flat  plate;  pressure  distributions  on  the  NACA  00 12 
airfoil  in  steady  flow;  and  steady  boundary  layers  at  zero  angle  of  attack 
[4,5].  In  1979,  the  unsteady  flow  field  was  studied  by  measuring  the 
excitation  velocities  induced  by  the  rotating  ellipse  [6].  In  1980  we  began 
the  first  extensive  effort  to  measure  the  unsteady  surface  pressures  and 
velocity  profiles  of  the  boundary  layer  on  the  airfoil  at  zero  angle  of  attack 
[6].  Continued  surface  pressure  measurements  and  boundary  layer  profiles  were 
taken  in  1981  with  the  angle  of  attack  varying  up  to  15  degrees  [3].  Wake 
measurements,  two-dimensionality  tests  and  numerical  boundary  layer 
calculations  were  performed  in  1982  and  continued  through  1983  (3).  Also 
beginning  in  1983  were  studies  of  the  surface  pressures  in  separated  flow  [7]. 
Wall  shear  measurements  began  in  1984  with  the  initial  design  of  the  hot  wire 
shear  probe  and  the  hot  film  tests. 

There  were  two  major  accomplishments  made  during  the  last  year:  first,  we 
were  successful  in  constructing  and  testing  a  hot  wire  type  surface  mounted 
shear  probe;  and  second,  we  were  able  to  use  a  traveling  wave  excitation  to 
generate  the  unsteady  flow. 

The  work  done  by  P.  ?.  Lorber  [3!  for  his  Ph.D.  thesis  included 
extensive  analysis  of  the  surface  pressures  and  boundary  layers  on  an  NACA 
0012  airfoil,  but  he  could  only  estimate  the  wall  shear  from  the  velocity 
profiles  of  the  boundary  layer.  To  vain  a  better  understanding  of  the  wall 
shear  associated  with  the  NACA  00 1 2  airfoil,  a  wall  shear  probe  was 
constructed.  The  probe  is  shown  in  figure  3  and  the  details  of  the  orobe 
construction  can  be  found  in  K.  J.  ? little' s  3.  M.  thesis  [1-.  This  probe, 
having  the  high  frequency  response  characteristics  of  a  hot  wire  type  device, 
enabled  us  to  measure  the  wall  shear  directly  for  steady  or  ’unsteady  flows. 
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The  other  major  accomplishment  was  creating  a  different  type  of  unsteady 
perturbation  and  measuring  the  induced  aerodynamic  behavior  from  it.  As  noted 
the  initial  the  unsteady  perturbation  was  created  using  a  rotating  elliptical 
cylinder  located  below  and  behind  the  airfoil  (see  figures  i  and  2).  The 
result  of  the  airfoil/ellipse  combination  lead  to  an  essentially  constant 
phase  type  excitation.  This  year  we  were  able  to  generate  a  traveling  wave  by 
means  of  the  gust  generators.  The  gust  generators  are  two  parallel  airfoils 
(see  figures  4  and  5)  which  create  a  periodic  circulation.  If  the 
circulations  of  the  two  generators  are  in  phase,  the  effective  perturbation  is 
in  the  longitudinal,  or  streamwise  direction.  If  the  circulations  are  out  of 
phase,  the  effective  perturbation  is  in  the  lateral,  or  spanwise  direction. 
The  details  of  the  gust  generators  performance  can  be  found  in  reference  8. 

TECHNICAL  ACTIVITIES 

SHEAR  PROBE  RESEARCH 

An  extensive  amount  of  wall  shear  data,  both  steady  and  unsteady  was 
recorded  during  the  last  year.  What  follows  is  a  summary  of  the  significant 
results  obtained  by  K.  J.  Flittie  and  reported  in  his  S.  M.  thesis  [1]. 

The  hot  wire  shear  probe  is  based  or.  a  probe  designed  by  J.  Cousteix  ’ 9 ] 
of  ONERA/CERT.  Essentially,  the  probe  is  a  hot  wire  located  flush  with  the 
surface  (see  figure  3)-  There  is  a  cavity  etched  into  the  surface  beneath  the 
sensing  wire  in  order  to  prevent  contact,  and  consequently  thermal  conduction, 
with  the  surface.  Thus  the  probe  retains  a  high  dynamic  response,  which  makes 
it  ideal  to  use  in  measuring  unsteady  flows. 

To  test  the  integrity  of  the  shear  probe,  the  turbulent  wall  shear  was 
measured  on  a  flat  plate  with  a  zero  pressure  gradient.  The  results  shown  in 
figure  6  demonstrate  that  the  measurements  from  the  probe  were  repeatable  to 
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within  101  of  the  predicted  empirical  results  using  boundary  layer  profiles. 
It  is  worth  noting  that  in  each  case,  the  measured  shear  stress  was  higher 
than  the  predicted  value.  It  is  believed  that  this  is  attributed  to  surface 
roughness  rather  than  a  property  of  the  probe,  as  will  be  discussed  below. 

A.  CONSTANT  PHASE  EXCITATION  WALL  SHEAR  RESULTS 

The  wall  shear  was  measured  for  steady  and  unsteady  flows  at  the  701,  851 
and  941  chord  locations  for  zero  and  ten  degrees  angle  of  attack  and  Reynolds 
number  of  700,000.  The  unsteady  cases  studied  consisted  of  four  different 
reduced  frequencies,  0.5,  1.0,  2.0,  and  6.4.  Unsteady  boundary  layer  profiles 
were  also  taken  for  the  above  cases.  There  were  four  significant  results 
found  during  this  program: 

(1)  The  time  mean  velocity  profiles,  integral  parameters  and  wall  shear 
stresses  were  the  same  as  in  the  equivalent  steady  cases  studied. 

(2)  A  quasi-steady  Clauser  method  was  developed  and  found  to  yield  time  mean 
c  's  to  within  151  of  the  experimental  cf's  and  phase  angles  within  50  degrees 
of  the  experimental  phase  angles  for  reduced  frequencies  of  less  than  2. 

(3)  For  unsteady  flow  near  the  wall,  a  universal  velocity  distribution  does 
not  appear  to  exist. 

(4)  The  oscillating  adverse  pressure  gradient  has  significant  influence  on 
the  viscous  sublayer. 

Figures  7  through  10  show  the  time  mean  u+  versus  y+  curves  at  zero  angle 
of  attack  and  851  chord.  Mete  hew  all  the  curves  are  very  nearly  identical. 
If  we  compare  these  curves  tc  the  steady  case,  figure  1  * ,  at  zero  angle  of 
attack  and  851  chord,  we  find  no  discernible  difference.  This  same  behavior 
can  be  seen  in  the  velocity  profiles  and  wall  shear. 
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The  method  proposed  by  F.  H.  Clauser  [10]  assumes  that  a  universal 
velocity  profile  exists  near  the  wall.  The  method  was  developed  for  steady 
turbulent  flow,  but  we  used  it  in  a  quasi-steady  analysis.  The  idea  is  to 
apply  the  Clauser  technique  at  various  points  over  a  cycle.  This  assumes  that 
the  logarithmic  overlap  layer  is  valid  at  each  point  in  the  cycle.  It  also 
assumes  that  the  friction  velocity,  u*,  is  the  correct  scaling  parameter  for 
the  inner  region  of  the  boundary  layer  and  that  the  velocity  phase  angle 
remains  constant  in  the  logarimthic  region. 

The  results  from  the  Clauser  analysis,  shown  in  figures  12-15,  are 
encouraging.  The  time  mean  cf's  from  the  Clauser  technique  were  within  I5t  of 
the  predicted  the  experimental  cf's  up  to  a  reduced  frequency  of  approximately 
2.  The  Clauser  technique  also  predicted  the  phase  angles  to  within  50  degrees 
of  the  experimental  results. 

It  appears  that  of  the  four  regions  identified  with  the  turbulent 
boundary  layer,  the  viscous  sublayer  is  the  most  effected  by  an  unsteady  flow. 
The  velocity  phase  shift,  for  example,  has  been  found  to  be  a  function  of 
reduced  frequency.  The  general  shapes  at  reduced  frequencies  of  above  and 
below  2  are  shown  in  figure  16  and  17.  It  is  important  to  note  how  phase 
shift  changes  rapidly  at  low  reduced  frequencies  as  the  wall  is  approached. 
The  hope  for  a  universal  velocity  distribution  is  not  substantiated  since  the 
assumption  for  such  a  principle  a  constant  phase  shift  as  the  wall  is 
approached,  fails. 

To  analyze  the  qualitative  behavior  of  the  phase  shift  with  reduced 
frequency,  consider  a  small  cer~ urbation  expansion.  The  ensemble  average 


momentum  equation  is 


(1) 


3<u> 

IF- 


+  <u> 


3<u> 

3x 


+  <v> 


3<u> 

3y 


1_  3<p>  3^<u>  3<u’u,>  3<u'u' > 

P  3x  V  3y2  3x  3y 


Since  we  are  considering'  flow  very  near  the  wall,  we  neglect  the  turbulence 
and  convection  terms  so  1  simpifies  to* 
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Where  Ue(dUe/dx)+dUe/dt  is  from  the  unsteady  Bernoulli  equation.  We  now 
expand  the  u  components  of  velocity  in  a  small  pertubation  parameter,  (e). 
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Substitution  3  and  4  into  2  and  collecting  powers  of  !  ■-  )  yields 
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Consider  the  real  part  of  the  first  order  term 
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The  first  term  on  the  right  hand  side  is  due  to  the  unsteady  external  intertia 
term.  The  second  term  on  the  right  hand  side  is  due  to  the  unsteady  pressure 
gradient.  For  small  reduced  frequencies,  the  first  term  (the  inertial  effect 
of  the  unsteady  external  velocity)  is  negligible  compared  to  the  unsteady 
pressure  gradient.  Clearly,  if 
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the  'unsteady  pressure  gradient  is  the  dominant  factor  in  determining  the  phase 
shift  and  amplitude  of  the  velocity  near  the  wall.  As  the  reduced  frequency 
is  increased,  and  if  the  unsteady  external  velocity  increases,  the  inertial 
term  becomes  significant  and  effects  both  the  phase  shift  and  amplitude. 


B.  TRAVELING  WAVE  EXCITATION  WALL  SHEAR  RESULTS 

During  the  traveling  wave  testing,  the  wall  shear  stress  was  measured  at 
the  same  test  Reynolds  number,  700,000,  angles  of  attack  and  chord  locations, 
70$,  85 Jf  and  94J.  Although  we  were  able  to  generate  both  a  longitudinal  and 

lateral  perturbation,  the  reduced  frequency  was  limited  to  1.5  due  to  the 
power  limits  to  drive  of  the  gust  generators  [8],  The  significant  results 
from  this  test  sequence  were  as  follows: 

(1)  As  before,  the  time  mean  velocity  profiles  and  wall  shear  values  were 
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the  same  as  the  equivalent  steady  cases  studied. 


(2)  The  lateral  perturbation  at  lower  reduced  frequencies  appears  to  create 
locally  high  pressure  gradients  which  accelerate  and  decelerate  the  flow 
resulting  in  highly  inflected  unsteady  velocity  profiles  and  an  oscillating 
transition  point  on  the  aft  portion  of  the  airfoil. 

(3)  The  quasi-steady  Clauser  method  showed  a  qualitative  agreement  to  the 
experimental  waveform  results  over  a  cycle.  Quantitatively,  as  before,  the 
time  mean  values  of  cf  were  within  15 t  of  the  experimental  values  and  phase 
shifts  to  within  50  degrees. 

The  time  mean  curves  of  0+  versus  y+,  figures  18-23,  demonstrate  how  the 
time  average  parameters  correspond  to  the  values  of  the  steady  cases. 
Comparing  figures  1 8-23  to  the  steady  curves  shown  in  figures  11,  we  find  no 
appreciable  differences.  Recall,  this  was  the  same  result  found  for  the 
constant  phase  excitation,  thus  we  conclude,  this  characteristic  is 
independent  of  the  type  of  excitation. 

Although  certain  characteristics  of  the  unsteady  boundary  layer  flows 
were  independent  of  the  type  of  excitation,  (that  is,  constant  phase  versus 
traveling  wave)  this  was  net  true  for  every  case.  For  example,  at  lower 
reduced  frequencies  in  the  lateral  mode,  there  was  a  strong  indication  of 

locally  high  pressure  gradients  and/or  laminar  to  turbulent  transition. 

Figures  24  and  25  show  the  boundary  layer  characteristics  for  the  70S  chord 

position  with  a  lateral  mode  excitation  and  reduced  frequency  of  0.5.  Notice 

that  the  skin  friction  coefficient  curve,  figure  25,  is  step-like  and 
considerably  smoother  between  phase  angles  of  90  and  240  degrees.  This  is 
very  different  than  either  the  longitudinal  case  (figure  26)  or  the  constant 


phase  case  (figure  27). 


Assuming  that  a  quasi-steady  analysis  is  valid  at  this  low  reduced 
frequency,  we  can  examine  the  flow  field  at  any  point  in  the  cycle  as  if  it 
were  steady.  A  steady  turbulent  boundary  layer  velocity  profile  becomes  more 
highly  inflected  with  increasing  adverse  pressure  gradient,  the  shape  factor 
increases  and  the  cf  decreases.  This  behavior  is  seen  in  figures  24  and  25 
and  it  appears  that  the  velocity  profile  oscillates  between  a  fully  developed 
and  highly  inflected  profile  on  the  aft  section  of  the  airfoil.  In  other 
words,  an  oscillating  transition  region  and  locally  high  adverse  pressure 
gradients  causing  highly  inflected  velocity  profiles  near  the  trailing  edge 
may  account  for  the  rapid  variations  seen  in  the  skin  friction  coefficient, 
shape  factor,  and  local  pressure  gradient  ensemble  average  plots. 

As  in  the  constant  phase  excitation  analysis,  a  quasi-steady  Clauser 
technique  was  used.  Figures  28-30  show  the  comparison  between  the 
experimental  data  at  the  70*  chord  and  the  quasi-steady  Clauser  method.  The 
ensemble  Clauser  wave  exhibits  a  qualitative  similarity  to  the  data,  with  best 
results  at  the  low  reduced  frequency.  The  time  mean  values  from  the  Clauser 
method  were  also  within  15 %  of  the  time  mean  experimental  values  and  the  phase 
shifts  from  the  Clauser  method  were  within  50  degrees  of  the  experimental 
phase  angles.  This  technique  proved  to  be  useful  for  both  types  of 
excitations,  yielding  common  results  for  each  case. 


UNSTEADY  PRESSURE  RESULTS  USING  A  TRAVELING  WAVE 
Four  years  ago,  for  his  3.  M.  thesis,  P.  E.  Lorber  [6]  used  -unsteady  thin 
airfoil  theory  to  predict  surface  pressures.  3y  measuring  the  upwash  induced 
by  the  unsteady  flew  along  ‘he  chord  of  the  airfoil.  Lorber  calculated  the 


surface  pressures  using  a  method  developed  by  Theodcrsen.  This  method  assumes 


a  constant  phase  excitation,  which  is  the  type  of  flow  field  generated  by  the 
airfoil/ellipse  configuration.  By  mearsuring  the  surface  pressures  directly, 
he  was  able  to  compare  his  theoretical  results  to  the  experimental  results. 

By  employing  a  traveling  wave  type  excitation,  Theodorsen's  analysis  can 
no  longer  be  used  to  describe  resulting  surface  pressures.  There  are  two 
theories  which  apply  to  a  convected  gust  (or  traveling  wave).  One  is  a 
variation  of  Theodorsen's  method  in  which  the  convected  gust  is  accounted  for. 
This  method  was  developed  by  Von  Karman  and  Sears  [  1 1  ] .  The  other  analysis , 
developed  by  Horlock  [12]  applies  to  longitudinal  gusts  over  a  symetric 
airfoil  at  small  angles  of  attack.  The  results  of  both  of  these  analysis  are 
described  in  S.  W.  Linton's  S.  M.  thesis  [2]  and  will  be  discussed  here. 

A.  LATERAL  GUST  RESULTS 

For  the  lateral,  or  vertical,  gust  cases,  S.  W.  Linton  [2]  applied  an 
extended  Sears  theory.  The  extended  theory  accounted  for  variation  of  the 
pressure  amplitude  over  the  airfoil  chord.  Using  approximations  this  theory 
to  model  the  unsteady  flow,  Linton  was  able  to  estimate  the  pressure 
distribution  on  the  airfoil.  By  comparing  the  estimated  results  to  the 
experimental  results,  he  found: 

(1)  The  magnitude  and  location  of  the  iifferer.ee  pressure  amplitude  lid  not 
model  accurately,  with  the  basic  theerv. 

f‘2)  The  extended  theory  protected  phase  variations,  tut  of  a  lover  nagr.etude 
and  more  gradual  than  the  iata. 

Figures  31-34  show  me  :  .r  it  tarmonic  :f  the  urmeady  difference 
pressure  coefficient,  the  'la.-.-i  -ars  theerv  prediction,  ar.i  the  extended 
Sears  theory  prediction.  A*  *-•>  .  -  moused  frequency  of  ' .  0 ,  teth  theories 
offer  practical  results,  t  .*  -.1  --'..'•m  frecuer.cy  increases,  both  theories 
breakdown.  The  classical  me: rv  :  not  seem  *0  predict  any  oscillatory 

behavior  above  a  reduced  frecu-t  -•••  *  r  t  *  .2* .  Alt  rough  the  extended  theory 
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does  predict  the  oscillations,  it  fails  near  the  leading  edge  and  trailing 
edge,  particularly  at  the  higher  reduced  frequencies  (see  figure  31).  This  may 
be  due  to  the  fact  that  phase  velocity  variations  were  not  modeled. 

Figures  35-38  show  the  phase  lags  over  the  chord.  According  to  the 
classical  theory,  the  phase  lag  should  be  constant  over  the  chord,  but  this 
was  clearly  not  the  case.  The  extended  theory  is  shows  some  variation  of  the 
phase  over  the  chord  is  to  be  expected  but,  the  expected  behavior  compares 
badly  with  the  measured  behavior  of  the  phase  lag.  These  results  seem  to 
indicate  that  the  pressure  waves  are  not  moving  along  the  airfoil  at  constant 
phase,  as  predicted  by  the  classical  theory. 

B.  LONGITUDINAL  GUST  RESULTS 

To  analyze  the  longitudinal  gust  cases,  theory  by  Horlock  [12]  was  used. 
This  is  a  higher  order  effect  chan  the  vertical  gust  excitations.  In  this 
case,  the  lift  is  proportional  to  the  angle  of  attack  times  the  longitudinal 
gust  amplitude.  We  see  how  small  this  effect  is  by  comparing  figures  39  and 
40.  Figure  39  is  the  full  difference  pressure  coefficient  due  to  a 
longitudinal  gust  at  a  reduced  frequency  of  0.5.  The  unsteadiness  is  very 
difficult  see,  but  by  subtracting  out  the  mean,  as  in  figure  40,  the 
unsteadiness  is  clear.  Because  this  effect  was  smaller  than  indicated  by  well 
to  the  theory  of  Horlock. 

T  E  0  H  N  1  0  A  L  3  U  M  M  A  R  Y 

Research  concluded  last  year  ur.ier  contract  ?!*Q20-79-026  with  further 
examination  of  the  turbulent  boundary  layer  and  surface  pressures  on  an  N’ACA 
0012  airfoil  subjected  to  an  unsteady  flow.  There  were  two  types  of  waveform 
used  in  generating  the  unsteadiness.  The  first  was  a  constant  phase  type 
excitation.  This  type  of  waveform  was  used  for  most  of  the  'unsteady  testing. 


The  second  type  of  waveform  was  a  traveling  wave  type  excitation.  Boundary 
layer  profiles,  wall  shear  measurements  and  surface  pressure  measurements  were 
taken  on  the  airfoil  when  subjected  to  the  traveling  gust  type  excitation. 
The  wall  sheer  measurements  were  successfully  made  at  the  70t,  85*  and 

airfoil  chord  locations  for  steady  and  unsteady  flows.  Both  travel  wave 
and  constant  phase  excitations  were  used  to  generate  the  unsteady  flow  fields. 
Boundary  layer  profiles  were  also  taken  at  these  locations  for  the  purpose  of 
comparing  the  measured  wall  shear  to  the  wall  shear  predicted  by  a  quasi¬ 
steady  Clauser  method.  Results  were  obtained  for  a  Reynolds  number  of 
700,000,  and  angles  of  attack  of  0  and  10  degrees.  At  0  degrees  angle  of 
attack,  the  experimental  results  were  within  15*  of  the  Clauser  predictions. 
At  10  degrees  angle  of  attack,  the  Clauser  method  failed  to  adequately  predict 
the  shear. 

Two  types  of  llnerarized  potential  thin  airfoil  theory  were  used  in  an 
attempt  to  model  the  traveling  wave  type  excitation.  For  lateral  perturbation 
cases,  the  classical  Sears'  theory  and  an  extended  Sears'  theory  were  used. 

For  our  excitation,  which  was  far  from  that  assumed  by  Sears,  the  classical 
theory  was  unable  to  predict  phase  or  amplitude  variations  over  the  airfoil 
chord.  The  extended  Sears'  theory,  which  allowed  for  the  amplitude  variation 
cf  the  wave  as  it  convects,  was  able  to  indicated  both  phase  and  amplitude 
variations,  but  results  ranged  from  adequate  at  lew  reduced  frequencies  '0.5) 
to  poor  at  the  higher  reduced  frequencies  ',.25-,.5).  For  longitudinal 
perturbations,  a  method  developed  by  Fori ock  was  osed.  However,  the  ’unsteady 
lift  developed  is  of  the  order  of  *r.e  angle  of  attack  multiplied  by  the 
magnitude  of  the  perturbation,  which  ma^es  this  a  small  effect. 

Consequently,  our  attempt  to  model  the  1  or.gitudir.al  gust  failed. 
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FIGURE  6  EXPERIMENTAL  AND  IHEORETICAL  FLAT  PLATE  KESL'LTS  FOR 
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